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Abstract. Adsorption of polymers to surfaces is crucial for understanding many 
fundamental processes in nature. Recent experimental studies indicate that the 
adsorption dynamics is dominated by non-equilibrium effects. We investigate the 
adsorption of a single polymer of length iV to a planar solid surface in the absence of 
hydrodynamic interactions. We find that for weak adsorption energies the adsorption 
time scales - iV(i+2'')/(i+''), where V is the Flory exponent for the polymer. We 
argue that in this regime the single chain adsorption is closely related to a field-driven 
polymer translocation through narrow pores. Surprisingly, for high adsorption energies 
the adsorption time becomes longer, as it scales ~ 7V'^^+''', which is explained by strong 
stretching of the unadsorbed part of the polymer close to the adsorbing surface. These 
two dynamic regimes are separated by an energy scale that is characterised by non- 
equilibrium contributions during the adsorption process. 
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Polymer adsorption is a fundamental phenomenon that controls many natural 
processes [1]. The adsorption of the polymeric molecules to different surfaces and 
interfaces is important for adhesion, colloidal stabilisation, development of composite 
materials and coatings, for cell adhesion and communication, and for protein-DNA 
interactions [1, 2]. The importance of polymer adsorption has motivated extensive 
experimental and theoretical investigations to understand the underlying mechanisms. 
As a result, the equilibrium properties of adsorbed polymers are now well-understood 
[3]. However, many experimental studies [4], supported by theoretical ones [5-7] 
indicate that non-equilibrium behaviour is increasingly important in polymer adsorption 
dynamics. 

One key parameter in polymeric adsorption is the height of the free energy barrier 
that monomers have to overcome in order to bind to the surface. If the barrier is high, 
one commonly calls the adsorption process chemisorption, while in the absence of a 
significant barrier, it is called physisorption. A further characterization of physisorption 
involves the strength of the binding interaction between each monomer and the surface. 
If this interaction is on the order of kT, the process is called weak physisorption, while 
one speaks of strong physisorption in the case of interactions of about 10 kT or more, 
as for instance typically encountered for hydrogen bonding. 

In chemisorption, the high barrier faced by monomers attaching to the surface 
slows down the adsorption process; this allows the adsorbed part of the polymer chain 
to partially relax in effectively equilibrium conformations, giving rise to formation of 
large loops via the accelerated zipping mechanism [7,8]. The absence of a significant 
barrier makes non-equilibrium effects even more important in physisorption [4,7]. It is 
not clear what mechanisms drive the polymer adsorption away from equilibrium in this 
regime [7]. One of the possible contributions is the interaction between neighbouring 
polymer molecules that can significantly slow down the overall dynamics. This source 
of deviation from equilibrium is commonly eliminated by considering the adsorption 
dynamics of single polymers [7,9-12]. 

The adsorption of single macromolecules for weak polymer-surface interactions has 
been investigated by a combination of analytical and computational methods [7,9-14]. 
Monte Carlo simulations with the bond fluctuation model revealed significant deviations 
from equilibrium dynamics [9, 10]. The adsorption time was reported to scale as 
~ fQj. self-avoiding polymers, while the exponent is equal to 1.50 ± 0.04 

when the excluded volume interactions are neglected. Computer simulations and an 
approximate theory were also used to investigate irreversible adsorption of tethered 
chains [13,14]. These investigations assumed that the polymer molecule during the 
adsorption has three parts: a segment of already bound monomers, a stretched linear 
part ("stem") and a remaining part ("flower") which is not affected by the force of 
adsorption. This theoretical model yields an adsorption time scaling as ~ A^" with 
a — 1 -|- 1/ pa 1.59. Here is the Flory exponent for the polymer, and i> ~ 0.588 
in three dimensions. Simultaneously, in the Monte Carlo simulations a smaller value 
of a, namely ^ 1.51, has been observed [14], but it was argued that finite-size effects 
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were responsible for this discrepancy. The stem-flower model was originally proposed 
by Brochard-Wyart [15] for polymer chains under strong flows (under constant and 
very large flow velocity). It has a clear physical picture that allows one to obtain 
specific predictions for the dynamical properties. However, the growth velocity of the 
adsorbed polymer has been shown to be not large [13,14], not constant, and in time 
it even decays to zero; hence the validity of the stem-flower model to adsorption in 
all situations is questionable. Thus, despite many attempts, mechanisms of the single- 
polymer binding to the surfaces are still not well-understood. In this paper we present 
theoretical arguments supported by simulation data that clarify several non-equilibrium 
features of single-polymer adsorption. 




Figure 1. Schematic view of the single polymer adsorption to a planar surface. 
Unfilled circles correspond to binding sites available for adsorption by monomers. Red 
filled circles describe the surface sites already occupied by the polymer. Black filled 
circles represent monomers that are not adsorbed on the surface yet. 

The polymer molecule during binding to the surface can be viewed as consisting 
of two segments: the adsorbed monomers and the block of free monomers not on the 
surface. Theoretical studies argue that adsorption (for weak interactions) can be viewed 
as a sequential zipping process [13,14] in which the size of the adsorbed block increases 
by one monomer a time. This sequential mechanism suggests that the single-polymer 
adsorption process is closely related to field-driven polymer translocation (detailed 
later in the text), which has been intensively studied in recent years [16,22]. During 
field-driven translocation, the polymer molecule moves through a pore sequentially, 
decreasing the number of monomers on the cis side of the pore and sequentially 
increasing the number of monomers on the trans side of the pore. Theoretical studies 
of translocation based on the microscopic dynamics of the polymer [16,17] showed that 
memory effects are crucial for understanding this process. The memory effects appear 
due to the finite time to dissipate away or replenish the local enhancement in the density 
of monomers at the pore. From these works on translocation, it is reasonable to expect 
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that the same memory effects in the polymer should also play a role in the adsorption 
of single polymers to a surface. 

Consider a single polymer molecule that is near a solid surface, which starts to 
adsorb to the surface as shown in Fig. [H We assume that there are uniformly distributed 
binding sites on the surface, that the energy of adsorption per site is equal to e, and that 
the distance between binding sites is the same as the size of each monomer. We then 
use a Monte Carlo based FCC-lattice polymer code in three dimensions for self-avoiding 
polymers, with the rigid flat surface placed at z = 0, and study the adsorption dynamics 
for a variety of polymer lengths and for different strengths of the adsorption energy, 
in the absence of hydrodynamic interactions. In this polymer model the individual 
monomers perform both reptation and "sideways" movements [18] with each kind of 
movement attempted with frequency unity, which provides us with the definition of time. 
This model has been used before to simulate the diffusion and exchange of polymers 
in an equilibrated layer of adsorbed polymers [19]. Recently, we have used this code 
extensively to study polymer translocation under a variety of circumstances [16,17,20]. 

The process of adsorption involves a change in the free energy AF for the polymer: 
the polymer loses energy due to the attractive interaction between the surface and 
the monomers, and loses entropy, as adsorption makes the polymer collapse into two 
dimensions from three. The value of the adsorption energy e dictates the sign of AF, 
and thereby determines the fate of a partially adsorbed polymer. At high adsorption 
energies, the polymer will nearly completely adsorb since adsorption is energetically 
favoured, while at low adsorption energies, it will desorb since desorption is entropically 
favoured. In between, there is a critical value e* at which a partially adsorbed polymer 
will both adsorb and desorb with equal probability, for which the entropy gain for 
desorption is exactly compensated by the energetic gain for adsorption. 

Consider a partially adsorbed polymer of length A^, with m monomers (counting 
from one end) completely adsorbed on the surface while the rest (A^ — m) monomers 
are moving freely (off the surface). If one assumes that the adsorbed part of 
the polymer takes the form of a self-avoiding random walk on the two-dimensional 
adsorbing plane, then the partition function of this polymer is given by Z = 

[Aa/i^mT^'"" -1] [As/if ""^(A^ - m)T^''''-i]. Here, 7^^^) = 49/32 and 7f^^ ^ 1.16 are 
two universal scaling exponents [21], and A2, A^, ^2 and /is are model-dependent 

, ,,,1,1 dAFie) ^ . , dAFie) 
quantities. For this partially adsorbed polymer — can be estimated as 



dm dm 
dAFie) 

—e + kBT\n(fis/fi2) + ^/N corrections. Equating — to zero then yields the critical 

am 

adsorption energy e* ~ fc^T In(/i3//i2) in the limit of A^ ^ 00. Note that this expression 

is only an estimate, since in the adsorbed state not all monomers of the polymer adhere 

to the surface; moreover, as has been demonstrated in Ref. [13], the adsorbed part of the 

polymer takes a very compact conformation — much more compact than a self-avoiding 

dAFie*) 

walk in two dimensions. Nevertheless, — = shows that e* is in the order of 

dm 

knT. 
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For our model we determine e* in the following manner. We start with a polymer of 
length N with N/2 monomers from one end constrained to the surface (i.e., constrained 
to 2; = 1) without an adsorption energy, while the remaining N/2 monomers are free, 
and equilibrate the polymer under this constraint (the free N/2 monomers encounter the 
surface only as a planar obstacle). At time t = 0, an adsorption energy e is introduced, 
and simultaneously the constraint is lifted. We repeat this exercise for polymer lengths 
ranging from = 100 to = 800, while tuning the suppression of the desorption rate 
by a factor of f = exp[— ^/(/ebT)], until on average the adsorbed part of the polymer 
neither grows nor shrinks. The results for the critical values v* for several polymer 
lengths are summarised in Table [H from this Table we conclude that v* ~ 0.34 ±0.01 for 
our model, and thus that e* /{ksT) = 1.08 ±0.03. Since we use e > 2kBT, our polymers 
always adsorb, and any reference to high or low adsorption energies will henceforth refer 
to e > e*. 



N 


V* 


100 


0.405 


120 


0.395 


140 


0.392 


160 


0.38 


200 


0.365 


400 


0.36 


800 


0.342 



Table 1. The critical desorption rate v* — exp[~e* / (kBT)] as a function of polymer 
length N. 

The specific manner in which we simulate surface adsorption is as follows. We take 
a polymer of length (A^ ± n*) with n* monomers from one end constrained to 2; = 1 
without an adsorption energy (a process we term "grafting" for later reference), and 
equilibrate the rest of the polymer in z > 0, i.e., during the equilibration process the A^ 
free monomers encounter the surface only as a planar obstacle. We index the monomers 
consecutively along the chain, starting with i = —n* for the grafted end. The free end 
is thus indexed by i = A^, and the last grafted monomer corresponds to i = 0. At 
t = we switch on the attractive interaction between the monomers and the surface, 
and simultaneously lift the constraint. The dynamics of the polymer for t > is then 
governed by, in addition to self-avoiding polymer dynamics, the fact that the ratio of 
probability of a monomer (including the grafted monomers) jumping from 2; = 1 to 
z = 2 and that of a monomer jumping from ^ = 2 to 2 = 1 is given by the Boltzmann 
ratio exp[—e/ (/c^T)]. Throughout this paper we choose n* = 30; since we use adsorption 
energies higher than 2kBT, this implies that the probability for the entire polymer to 
detach from the surface is practically zero. It should be noted also that the specific 
value of n* does not affect the adsorption dynamics as long as n* <^ A^. 
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Given this setup, on average we expect the monomers to be adsorbed on the surface 
in a sequential zipping manner: on average monomer rii (> 0) will be adsorbed on the 
surface {i.e., attain z = 1) for the first time earlier than monomer n2 > ni. For future 
reference, at any time t for any configuration of the adsorbing polymer we can identify 
the monomer with the highest index n{t), which has 2; = 1, to be called the "active 
monomer". This definition divides the entire polymer into two segments: (i) a part 
consisting of monomers i < n{t), largely adsorbed to the surface, and (ii) another part 
consisting of monomers i > n{t) that behave as a polymer of length [N — n{t)] tethered 
on the surface at the location of the active monomer. 

In light of n{t) defined in the above paragraph, it is important to note that our 
setup involving the initial grafting of the polymer — albeit simplified — does capture 
the adsorption dynamics in a real situation. In reality, a long polymer does not adsorb 
starting from one end; almost always it starts to adsorb somewhere in the middle. 
Imagine, for a polymer of length N, a situation when the monomer with index uq 
(somewhere in the middle of the polymer), is the first one to adsorb. By definition, 
this monomer immediately divides the polymer into two separate "sub-polymers" — 
of lengths no and {N — no) respectively. If £ > £*, then (on average) these two sub- 
polymers will start being adsorbed independently (in stating this, we disregard the steric 
interactions between them) from their common end: one from monomer uq towards 
monomer 1, and the other from monomer uq towards N. The adsorption dynamics for 
a polymer in a real situation — at least in the scaling sense, which is the main focus of 
this paper — is the same as that of the polymer in our setup (which starts to adsorb 
from one end). Our setup — similar to the existing ones [13, 14] — therefore, is purely 
a choice of convenience to study the adsorption dynamics in a real situation. 

Returning to our setup, if adsorption were a sequential zipping process for every 
single realisation, then the adsorption dynamics can be described solely by the active 
monomer index n{t) as a function of time, and the dynamics of adsorption can be 
mapped exactly on to that of field-driven translocation. More precisely, in polymer 
translocation driven by a potential difference AV across the pore, when a monomer 
crosses from the cis (trans) to the trans (cis) side, the length of the polymer segment 
on the cis side reduces (increases) by one monomer with an energy gain (penalty) of 
magnitude qAV, where q is the charge of one monomer. Similarly, (a) if the active 
monomer happens to detach from the surface (with an energetic penalty e) then the 
length of part (ii) of the polymer increases by roughly one monomer; (b) alternatively, 
if the index of the active monomer increases by one (with an energy gain of s), then 
the length of part (ii) of the polymer decreases by one monomer. In Ref. [16], based 
on memory effects in polymer dynamics, two of us showed that the total number of 
translocated monomers at time t increases as a power- law ~ t(i+^)/(i+2^) at weak fields; 
recently, this has been confirmed by a different polymer models [22,23]. This implies that 
if the adsorption process were a sequential zipping process for every single realisation 
for our setup, n{t) would scale ~ Based on this result — although in a real 

situation adsorption is a sequential zipping only on average, and not for every single 
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realisation — on average we expect n(t) to increase in time t also as ~ 

for our setup, i.e., the adsorption time scales as ~ ^ when the adsorption 

energies are not very high. We demonstrate this in the paragraphs below. 

It is important to note that during adsorption for our setup, a monomer with index 
[n{t) + ni] may get adsorbed with the surface before any of the in between monomers 
[with indices n{t) + 1, . . . ,n{t) + Ui — 1] do. For such an event, the adsorbed part of 
the polymer is said to form a "loop" of length ni between monomers with indices n{t) 
and [n{t) +ni]; in fact, it is precisely such "loop formations" that prevent adsorption — 
unlike translocation, for which the first passage of the monomers through the pore takes 
place strictly sequentially — from being a sequential zipping process for every single 
realisation. Consequently, the traditional way to follow the progress of adsorption for 
our setup is to track the average total number of adsorbed monomers s{t) at time t, 
so that s{Tad) ~ N would define the adsorption time Tad- However, since s(t) for any 
single realisation will saturate at a value ~ 0{N), care needs to be taken in measuring 
s{t), otherwise saturation effects might affect the numerical determination of the true 
exponent. In order to avoid saturation effects, we define tn as the average time, and s„, 
as the average number of adsorbed monomers when the n-th monomer attains z = 1 for 
the first time, with the condition that no monomer with index > n has ever attained 
z = 1. Since t„ is defined only till n = N, this method ensures that s„ never saturates. 
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Figure 2. Scaling of s„ as a function of n, for several values of N and e; all curves 
collapse on a single master curve, corresponding to the scaling s„ '~ n, represented by 
the solid black line. 

Indeed, we find that adsorption of the individual monomers is not a sequential 
process for every single realisation precisely because of the loop formations as discussed 
in the above paragraph; however, as shown in Fig. [21 s„ does scale linearly with n, 
confirming that adsorption on average is indeed a sequential zipping process. This 
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Figure 3. (a) Weak adsorption data for N = 1000 (from left to right) e = 2 (black), 
£ = 3 (blue), e — 4 (magenta) and e = 5 (orange): the data are progressively 
separated by a factor 2 along the x-axis for clarity. The original data are shown in 
the inset in the same colour scheme. The solid black line corresponds to an exponent 
(1 + + 2v) ~ 0.73. (b) Strong adsorption data for TV = 1000 (from left to right) 

e = 6 (red), e = 7 (blue), e = 8 (magenta), e — 9 (orange), e = 10 (brown) and e = 25 
(green): the data are progressively separated by a factor 1.5 along the x-axis for clarity. 
The original data are shown in the inset in the same colour scheme. The solid black 
line corresponds to an exponent 1/(1 + v) ~ 0.63. Data correspond to an average over 
~ 400, 000 realisations for each value of N and e. Right panel: Collapse of the data for 
s{tn) VS. tn for e = 3, and TV = 200 (red), N = 500 (blue) and iV = 1000 (magenta), 
(c) Comparison of adsorption speed at several energies and A^-values: N = 1000, £ = 2 
(red), N = 1000, £ = 3 (blue), N = 500, £ = 3 (cyan), N = 200, £ = 3 (magenta), and 
— 1000, £ = 10 (black). Note that adsorption is slower for large e than for small e. 

property ensures that the exponent we get for s„ as a function of t„ is the same as 
the one that one would get from tracking s{t) as a function of t. A remarkable feature 
of Fig. [2] is the collapse of all s„ vs. n curves on a single master curve: it shows that the 
proportion of monomers in the loops within the adsorbed part of the polymer, given by 
(n — s„), is independent of e, a feature that we will return to shortly. 

In Fig. [3](a) we present the data for weak interactions with the surface, i.e., when 
the adsorption energy is not too high {e < 5), for which we do obtain the exponent 
(1 + z/)/(l + 2z/), corresponding to Tad ~ A^(i+2'')/(i+''). Additionally, the data exhibit 
energy-dependence [see the inset, and also Fig. [3](c)], demonstrating that for e < 5 the 
higher adsorption energy also yields faster adsorption, like higher field means shorter 
(field-driven) translocation time at weak fields [16]. The situation changes for stronger 
interactions {e > 5): in Fig. [3](b), we register a slowly decreasing slope in the t„- 
Sn log-log plot with increasing adhesion energy; eventually for the virtually irreversible 
adhesion process e = 25, we recover an exponent 1/(1 + 1/), i.e., Tad ~ A^^^'^, in agreement 
with Refs. [9, 10, 13, 14]. For these values of e, Tad is independent of e. 

The surprising aspect of two different scaling regimes for Sn vs. tn, as shown in Figs. 
[3]^a) and (b) is that for long enough polymers adsorption is faster for low £- values than for 
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high £- values. This is demonstrated in Fig. El^c). In fact, Fig. [3](c) leaves one wondering 
whether the slowdown of adsorption is due to phenomena at the adsorbing surface that 
are different for strong adsorption energies than for weak adsorption energies. E.g., for 
high e-values it has been shown in Ref. [13] that the adsorbed part of the polymer takes 
a very compact form. Based on this result of Ref. [13], it may be argued that since for 
high £- values the individual monomers are essentially irreversibly adsorbed, the polymer 
needs to form systematically bigger loops to access available surface sites for adsorption, 
a phenomenology that is absent for low e-values; and as a result adsorption is slower 
for high e-values than for low er-values. Such a possibility is however ruled out by the 
collapse of the data over a very wide range of e values in Fig. [21 it shows that on average 
the fraction of monomers in the loops [given by {n — Sn)] has no dependence on e; i.e., 
steric hindrances due to the adsorbed part of the polymer does not cause the slowdown 
of adsorption at high e-values. Instead, as explained below quantitatively, the physics 
of the slowdown of adsorption at larger e-values is explained by the lack of availability 
of not-yet-adsorbed monomers near the surface. 

For high adsorption energies the monomers that were close to the surface at t = 
initially get quickly and effectively irreversibly adsorbed, while the monomers that are 
far away from the surface at t = cannot respond to this fast change of the polymer's 
configuration near the surface. As a result, during the adsorption process, the polymer 
adopts a stretched configuration close to the surface, while far away from the surface 
the polymer remains largely in its t = coiled shape: this is the stem-flower picture of 
Brochard-Wyart [15]. [It is precisely this stem-fiower shape that invalidates the physics 
behind the exponent {1 + + 2z/), seen at low adsorption energies. The number 

(1 + 2z/) in the denominator is derived from the Rouse exponent, and the number 
(1 + z/) in the numerator assumes that during adsorption the polymer's size scales as 
~ [N—n{t)Y; both fail at the stem (of the stem-fiower model), which is highly stretched.] 
In fact, the occurrence of u in the exponent 1/(1 + z/) at high adhesion energies stems 
from the polymer's size-scaling N'^ at t = 0, as we argue next. Let us denote, by z{t), 
the distance that the stem extends in real space from the surface at time t. The total 
number of monomers in the flower at time t — still largely in the same coil shape as at 
t = — is Q{t) ^ N — z{ty^'^ . In such a configuration, the flower would lose monomers 
through the stem to the surface, and the rate of loss of monomers is proportional to the 
gradient of monomeric density along the stem, ~ l/z{t). The solution of the differential 
equation Q{t) ~ ^/z{t) yields z{t) ~ t'^/i'^+'^). Since all the monomers that were present 
within a distance z{t) at t = — apart from the few within the stem at z{t) — 
are adsorbed by time t, the total number of adsorbed monomers at time t scales as 
s(t) ~ z{ty^'^ ~ Note that in this qualitative derivation there is no dependence 

on the adsorption energy [except that it needs to be high!], as observed in the inset of 
Fig. [3](b). When hydrodynamic interactions are included, following the physics of field- 
driven translocation [16] we expect the adsorption time to scale ~ iV^^^'^^/^^'') for not 
too high adhesion energies; however, presently we do not understand how co-operative 
motions of the monomers in the presence of hydrodynamics would affect the exponent 
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at high adhesion energies. 

An increase in the energy of adsorption allows one to cross-over from weak to strong 
regimes of physisorption, both characterised by different exponents. The remaining 
question is what determines the energy scale Ec that separates these two regimes. This 
can be understood if we return to AF ~ — e + ksT In [fi^/ ^12), wherein the first term 
lowers the free energy due to favourable adsorption to the attractive surface, while 
the second term increases the free energy because of entropy reduction by going from 
three dimensions to a more constrained two-dimensional surface. It is reasonable to 
suggest that two dynamic regimes of adsorption are separated when the free energy 
gain per monomer is comparable with thermal energy, i.e., |AF| ^ ksT. Our estimates 
for critical adsorption yield ;U3/yU2 ~ 3, which leads to Ec — 2. Our simulations show 
that Ec ~ 5, suggesting deviations from (equilibrium) free-energy concepts in (non- 
equilibrium) surface adhesion process. As shown in Fig. EJ^c), the adsorption proceeds 
faster for lower adsorption energies, with the most optimal adsorption speed close to Ec- 
Since these type of energies are typical for protein-DNA interactions [2] , one can suggest 
that this might be a mechanism by which biological adhesion processes are controlled. 

To conclude, using computer simulations and theoretical arguments we studied 
single polymer adsorption to solid surfaces in the absence of hydrodynamic interactions. 
Our analysis shows that the adhesion process is non-equilibrium. Details of the 
adsorption process depend on the strength of adsorption energies: for weak (polymer- 
surface) interactions the dynamics is determined by memory effects as in field-driven 
polymer translocation, while for strong interactions adsorption is limited by stretching 
of the unadsorbed part of the polymer. These two regimes are separated by the energy 
scale that is determined by a balance between favourable enthalpic and unfavourable 
entropic contributions due to adsorption of the monomers to the surface. It is argued 
that the adsorption process is most optimal at low interaction energies, and this might 
be the mechanism by which biological surface adhesion processes are controlled. 
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